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^ ■ Abstract 
O 

Selfconsistent Hartree-Fock-Bogoliubov (HFB) calculations have been per- 
' formed with the Gogny force for nuclei along several constant Z and constant 

N chains, with the purpose of extracting the macroscopic part of the binding 
^ ' energy using the Strutinsky prescription. The macroscopic energy obtained 

in this way is compared to current liquid drop formulas. The evolution of the 
^ ' single particle levels derived from the HFB calculations along the constant Z 



and constant N chains and the variations of the different kinds of nuclear radii 
are also analysed. Those radii are shown to follow isospin-dependent three 
parameter laws close to the phenomenological formulas which reproduce ex- 
perimental data. 
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I. INTRODUCTION 



Microscopic calculations using the effective nucleon-nucleon interaction proposed by 
Gogny have been shown to reproduce in very satisfactory way a variety of nuclear prop- 
erties over a wide range of proton and neutron numbers For instance, binding energies, 
separation energies of one or two nucleons and charge and neutron mean square radii usually 
agree quite well with known experimental data. In view of this, it seems appropriate to use 
theoretical results obtained in such a framework to derive systematics of nuclear properties 
that could be useful for making predictions concerning nuclear species or aspects of nuclear 
structure no yet known as neutron distribution radii or binding energies of superheavy nuclei. 

In the present work, simple formulas are used in order to represent the binding ener- 
gies and the charge, proton and neutron root-mean-square (r.m.s.) radii obtained in HFB 
calculations with the Gogny force, and they are compared with existing phenomenological 
formulas. Our aim is to determine simple but realistic enough parametrizations that would 
avoid performing time-consuming microscopic calculations when only rough estimates of the 
binding energies or radii are needed. Some preliminary results have already been published 
in P,^, for a smaller number of nuclei however. 

Section |I| presents a short overview of the theoretical models used in this work. Sec- 
tion |T| explains how the basis parameters employed in the self-consistent calculations have 



been determined. In Section the results concerning the single particle levels schemes of 
representative isotonic and isotopic chains are analysed. Section gives the method used 
to substract the neutron and proton shell and pairing corrections from the HFB energy. 
The "macroscopic Gogny energy" obtained in this way is then parametrized by means of 



a liquid-drop like formula In Section ^I|, systematics of the mass, charge, proton and 
neutron radii and of the ratios of proton to neutron radii are presented. Simple isospin- 
and mass number- dependent parametrizations of these quantities are determined, which ap- 
pear in excellent agreement with those derived using the relativistic mean field method 0. 
The formula obtained for the deformation-independent ratios of proton to neutron r.m.s. 
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radii is especially useful, since it should be valid for all even-even nuclei, either spherical 
or deformed. Conclusions and perspectives for further work are gathered at the end of the 
paper. 



II. THEORY 



The Gogny density-dependent effective nucleon-nucleon force is of the following form [|l| : 

2 (fi - r2f 

i=l 

+ i M^Ls (Vi - V2) X 5{fi - f2)(Vi - V2) ■ (ai + ^2) (2.1) 

,ri + f2. 



+ to (1 + xqP,) 6{n - r2) 



Pi- 
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which represents a central finite range interaction, a zero-range spin-orbit term and a zero- 
range density dependent interaction, respectively, to which one has to add the Coulomb 
interaction in the case of protons. The central interaction is made up of two distinct gaussians 
with ranges fii and ^2 whose values are given below. Pa and Pr denote the spin and isospin 
exchange operators respectively, and p represents the total density. 
We use the Gogny DIS interaction the parameters of which are given below: 

Wi = -1720.30 MeV W2 = 103.639 MeV 

Bi = 1300.00 MeV B2 = -163.483 MeV 

Hi = -1813.53 MeV H2 = 162.812 MeV 

Ml = 1397.60 MeV M2 = -223.934 MeV (2.2) 

/ii = 0.7 fm /i2 = 1.2 fm 

to = 1390.6 MeV fm^(^+^) xq = 1 

7= 1/3 Wls= 130 MeV fm^ 

When pairing correlations are neglected, the HFB approach reduces to the Hartree-Fock 
(HP) method which determines a self-consistent approximation of the nuclear mean-field. 



The corresponding ground state energy Ehf includes a contribution i?shcii from shell effects 
which can be evaluated by applying the Strutinsky smearing procedure [0] to the HF single- 
particle level distribution. The remaining part of Ehf can be considered as a macroscopic, 
liquid- drop like contribution E. Denoting by e^/ the HF proton or neutron single-particle 
energies, the Strutinsky shell correction energy -Esheii is- 

-Eshell = X! ^'^1' ~ ^ ' (^■^) 

V 

where E is the smoothed energy 

E = 2 ep{e)de. (2.4) 



With the Strutinsky prescription, the smoothed density p is given by 



,e — e 



i2 



Pie) = Pie') e 7 /(^) de' , (2.5) 



where A is the Fermi energy, / the Strutinsky (6*^^ order) curvature correction polynomial 
and 7 the width over which smoothing of the single particle level scheme is performed. The 
reliability of this procedure requires that -Esheii displays a plateau when drawn as a function 
of 7, separately for protons and neutrons. A crucial parameter in this respect is the number 
of single particle levels taken into account in the evaluation of the shell correction. This 
problem will be discussed in Section 0. 

Subtracting from the Hartree-Fock energy the shell correction energies of neutrons -E^heii 
and of protons -Efheib S^^s the following estimate of the macroscopic part of the binding 
energy: 

-E'macr = EhF - E^hell " -^fhell • (2-6) 

Another correction to the macroscopic part of the binding energy is the contribution 
-Epair of pairing correlations. This correction has been calculated for each nucleus by taking 
the difference between the energy Ehfb obtained in a full HFB calculation with the Gogny 
force and the energy Ehf computed using the simple HF method: 



-E-pair — Ehfb — Ehf ■ (2.7) 

The obtained r.m.s. radii of the neutron, proton, charge and mass density distributions 
are obtained in the microscopic HFB calculations in the usual way. In deformed nuclei, these 
different radii contain the effect of the quadrupole and other deformations However, the 
ratios of proton to neutron radii are almost deformation-independent since the proton and 
neutron density distributions are very similarly deformed. 

Starting from the HFB mean square radii (r^), equivalent spherical liquid-drop radii R 
can be defined through the expression: 



R=\hi^)- (2-8) 



which follows from the formula giving the mean square radius of a uniform spherical density 
distribution of radius R: 

{r') = . (2.9) 

These equivalent radii have then been fitted to a three parameter formula, with explicit 
isospin and A dependence 

R = roo(l + + ^)A^ = ro{A, I)A^ . (2.10) 

similar to the formula used in Ref. ||^ in the analysis of the radii obtained from relativistic 
mean field calculations. As the ratios of the proton and neutron r.m.s. radii are almost de- 
formation independent, they could be well approximated by the following formula depending 
only on the numbers of protons and neutrons: 

^=c(l + a^^ + A). (2.11) 

rn A A 

This last relation, together with the effect of the non-point like charge distribution of the 
proton evaluated with the approximate formula: 

{rl) = {rl) + OMim\ (2.12) 



can be useful for estimating the neutron radius of a nucleus when its r.m.s. charge radius 



rch = V {rch) is known: 



;(rg,)-0.64fm^ 



III. PARAMETERS OF THE CALCULATION 

The nuclei studied in the present work are represented by crosses and dots in the {N,Z) 
plot of Fig. I. These are the nuclei close to magic proton and neutron numbers and along 
the /3-stability line whose ground states are expected to be spherical flOl- They include Ca, 



Sr, Sn, Sm, Pb and Th isotopes, the A^=50, 82 and 126 isotone chains between the proton 
and neutron drip lines and a few /9-stable spherical nuclei between A = 38 and A = 218. 

In the microscopic HF and HFB calculations, the self-consistent equations have been 
solved in matrix form by expanding the single particle or quasiparticle states on finite bases 
made of spherical harmonic oscillator (HO) eigenf unctions. These bases depend on two 
parameters: the number No of major HO shells included in the bases and the oscillator 
parameter hu. 

In the present study, bases with No=12, 14, 16 or 18 major shells have been taken into 
account, depending on the nucleus under study, the criterium being that No is large enough 
to ensure convergence of the HFB energy (-Ehfb) within a few keV. For each nucleus and 
each value of No, the parameter hu has been chosen as the value Tiujmin that minimizes the 
HF energy -Erf- The values of humin obtained for Ca, Sn and Pb isotopes, and for the 
A^=126 isotones are plotted in Fig. ^ as a function of A for A''o=14, a value of No large 
enough for all these nuclei. The error bars indicate the ranges of huj for which the variation 
of Ehf does not exceed 100 keV. One can see that, with this number of HO shells, relatively 
large changes in hu do not significantly affect calculated values of Ehf- Consequently, 
approximate analytical interpolation formulas can be used. 
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A large scale investigation of the values found in spherical nuclei shows that an interpo- 
lation formula such as: 

hoomin = n(l + k ^ (3.1) 

can be adopted for all even-even nuclei, with n and k depending on the number No of shells 
included in the bases. In the present work, all the obtained values of humm could be fitted 
with the values n=64.05 MeV and fc=-0.46. A even more accurate interpolation formula, 
represented by the solid curve in Fig. |^, has been obtained for the specific set of nuclei for 
which No=14: is found to be a sufficiently large number: 

hu^in = (0.0002^2 - O.IA + 21.1) MeV (3.2) 

This formula gives a better agreement to calculated Tiojmin for the Ca up to the Pb isotopes 
and for the A^=126 isotones. 



IV. SINGLE PARTICLE LEVELS 

In order to better visualize the origin of shell effects in the nuclei considered in this 
work, the Hartree-Fock single-particle energies obtained with the Gogny force are plotted 
in Figs. § and ^ for neutrons (left hand sides) and protons (right hand sides). The 
three figures correspond to A^=50, A^=82 and A^=126 isotones, respectively, and the levels 
are drawn as functions of the proton number Z. They are labelled with the usual radial 
quantum numbers n, orbital angular momentum I and total angular momentum j. The 
numbers Ns between the levels indicate the number of nucleons necessary to fill all the 
levels located below. 

One can see that, in addition to the magic shells corresponding to Ns= 20, 28, 50, 82, 
126 and 184, well-marked subshells appear depending on the number of neutrons N . For 
instance, subshells with iV5=40 are clearly visible in the neutron and proton level schemes of 
the A^=50 and A^=82 isotones, and Ns=Q^ also appears as a subshell in the iV=126 isotones. 
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In Figs. 1^ to |TT] the neutron (left hand sides) and proton (right hand sides) single particle 
levels of the Ca, Sr, Sn, Sm, Pb and Th isotopes are plotted as functions of the neutron 
number A^. Subshells corresponding to A^5=40, 64 or 100 can also be observed in these 
figures besides the usual magic numbers A^5=50, 82, 126, 184. It is interesting to note that 
the magnitude of the proton and neutron A^5=50 shell gaps steadily decreases from ~ 6 
MeV down to ~ 4 MeV as the proton number Z increases from 20 to 90. This decrease does 
not occur in the case of the other magic numbers Ns=82 and A^5=126. 



V. MICROSCOPIC AND MACROSCOPIC PART OF BINDING ENERGIES 

As explained in Section the HF single-particle levels shown in Figs. |^ to |ri| have been 
used to separate HF binding energies into two components by means of Strutinsky's method: 
a macroscopic part behaving smoothly as a function of A and a shell correction reflecting the 
neutron and proton shells. The method of Strutinsky is a subtle technique which demands a 
very careful choice of the number of single-particle states above the Fermi level included in 
the smoothing procedure so as to obtain a plateau in the variations of the shell corrections 
with respect to the smoothing width 7 . In particular, including too many single particle 
states with positive energies strongly affects the determination of the shell effects associated 
with occupied orbitals. This is so because single particle states with positive energies are 
obtained in the present microscopic approach as a discrete set whose energy spectrum and 
shell structure can strongly depend on the choice of basis parameters. 

For this reason, the HF single-particle proton and neutron levels included in the Struti- 
nsky smoothing procedure have been chosen in the following way: For lighter nuclei, a cut 
off in the single particle energy has been introduced, having the usual value A + 2hujQ, with 
A the proton or neutron Fermi energy and huJo the average spacing between major shells. 
The number of single-particle levels taken into account in this way is of the order of 30, with 
maximum energies 5 MeV for neutrons and 15 MeV for protons above the Fermi level. In 
heavier nuclei, the number of levels included in the smoothing procedure was increased to 



8 



50, which corresponds to a cut-off in the single particle energies of the order of X + Ahuo, i.e. 
15 MeV for neutrons and 30MeV for protons. These numbers have been found to provide a 
satisfactory stabilization of shell corrections in all the nuclei studied in this work. 

In order to find a plateau in the variations of the shell correction with respect to the 
width 7 appearing in eq. (^2.51), 7 has been varied, taking as a scale the average shell spacing 



hujo = 4:0 A 3 MeV fTll]. In fact Strutinsky calculations generally obtain a plateau in the 



shell correction energy of both protons and neutrons for 7 in the vicinity of the traditional 



value 7 = 1.2huJo [|ri|]. Fig. |T2| and give examples of the variations with 'y/{fiujQ) of 
the shell corrections obtained in the present work. Similar behaviours are found in all the 
nuclei displayed in the {N,Z) plot of Fig. |l|. In Fig. [1^, the neutron (solid lines) and proton 



(dashed lines) shell corrections of twelve N=82 isotones ranging from Nd to Pb are plotted 
as functions of 'y/{huJo). In Fig. |13|, the different curves represent the shell corrections for 
neutrons (top) and protons (bottom) calculated in N=82 (left) and N=126 (right) isotones. 
One can see that, except for a few cases, as the proton shell corrections in ^^^Yh and ^^^Dy, 
reasonable plateaus are always found around the value 7 = 1.2 huQ. The results shown in 
the next figures are those obtained with the latter value of the averaging width 7. 

In the upper part of Fig. ffl the neutron and proton (dashed lines), and total (solid lines) 
shell corrections are drawn as functions of A for the Ca up to Th isotopes (left), A^=50, 82 
and 126 isotones (center) and /3-stable nuclei (right) included in this study. The lower part 
of this figure displays the corresponding pairing corrections calculated from eqs. ( |2.7|) . One 
observes that the variations of the shell corrections nicely reproduce the expected pattern, 
with minima located at magic proton and neutron numbers. Neutron shell corrections are 
seen to be almost independent of the nucleus neutron number A^, as they should, and proton 
shell corrections are found practically independent of Z. 

The macroscopic energies ( |2.6| ) obtained by subtracting the proton and neutron shell 
corrections from the Hartree-Fock energies Ehf for all the nuclei displayed in Fig. |I] are 
represented by lines in the three upper plots of Fig. |r^. A fit of these macroscopic energies 
using a liquid drop formula of the form proposed by Myers and Swi§tecki yields the 
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following expression: 



^2 '71 

^LD = [15.65(1 - 1.92/2)v4 - 18.92(1 - 2.1/2)^2/3 - 0.73-— + 1.99—] MeV . (5.1) 

with / = (A^ — Z')jA. The numerical parameters in this formula are in good agreement 
with those of Myers and Swicitecki (M-S) 0, except for the coefficient in the last term: 1.99 
instead of 1.21. This term represents the correction to the Coulomb energy of a charge liquid 
drop that accounts for the diffuseness of the proton density distribution. The correction 



appears to be larger in the case of formula (|5.1| ) - the difference with the equivalent M-S 
term reaches ~ 30 MeV in ^^^p,^ _^ which compensates for the fact that the Coulomb energy 
derived from the self-consistent results is larger than the M-S one. Let us also note that 
the asymmetry coefficients in the volume and surface terms of (|5.1| ) - the coefficients of 
/2 - are not equal, contrary to what is assumed in the M-S formula. On the other hand, 
the parameters appearing in (|5.1| ) are closer to those of the most recent best empirical fits 



- where the whole set of presently known nuclear masses and fission barriers is used - than 



to those of the various formulae derived from liquid drop models |]T2|. 

The fitted macroscopic energy ( |5.1| ) and the calculated values -Emacr are completely su- 
perimposed in the three upper plots of Fig. ITR In order to better appreciate the quality of 



the fit given by eq. (|0| ), the differences -Emacr — -Eld are displayed in the lower part of the 
figure, using an enlarged energy scale. One observes that these differences do not exceed 3 
MeV, which represents less than 0.5 % of the total Euf energy. 



VI. RADII 

In order to derive systematics for the radii of all the nuclei displayed in Fig. [l|, eq. ( p.8|) 
has been used to extract from the different - neutron, proton, mass and charge - Hartree- 
Fock-Bogoliubov r.m.s. radii, corresponding isospin-dependent radius constant Tq, Tq , Tq"* 
and r^. Namely, each kind of microscopic r.m.s. radius is multiplied by -y/|A~3. A fit of 
the obtained values with formula ( |2.10| ) yielded the following parametrizations: 
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= 1.17 (1 + 0.12J + 3.29/A) fm , 



(6.1) 



= 1.21 (1 - 0.14/ + 1.83/ A) fm , (6.2) 
= 1.19 (1 + 0.03/ + 2.70/ A) fm . (6.3) 



1.22(1 - 0.15/ + 2.32/A) fm, (6.4) 



„ch 



It is interesting to observe that the leading coefficient in ( p. 3D is almost equal to the usual 
value ro= 1.2 fm, although the / and A dependence included in the next terms may lead to 
significant deviations in the case of light or exotic nuclei. 

The above parametrizations appear very close to those derived in |^ by taking into 
account the Ca, Sr, Sn, Sm and Pb isotopes only, and to those obtained in the analysis of 
the r.m.s. radii calculated in the framework of the relativistic mean field theory for a set of 
nuclei similar to the one envisaged here P]. 

Using formula ( p.ll| ), the ratios of proton to neutron radii could be fitted with the 



parametrization: 

^ = 1.04(1 - 0.27/- 1.12M) . (6.5) 

This parametrization again appears very close to those derived in the studies of Refs. 
mentioned above. 

Let us point out that the parametrization given by formula ( |6.5|) , which has been es- 
tablished for the spherical nuclei shown in Fig. |^, has been found later on to reproduce the 
ratios of the HFB proton to neutron radii of all nuclei between the proton and neutron drip 
lines, spherical or deformed, in the range 20 < A < 318. This latter result follows from large 
scale HFB calculations with axial symmetry presently under way with the Gogny force. 

The results obtained in the present work are summarized in Figs. |1^ and 0. The upper 
part of Fig. displays the mass HFB radii (solid lines) at constant Z (left), constant 
(center) and for /^-stable (right) nuclei, together with the fits given by eq. ( |6.3|) . The lower 
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part of the figure sliows in a similar way the results obtained for charged radii. The fits 
are given by eq. ( |6.4| ), and experimental data, taken from Ref. [|1^], have been indicated by 
crosses. Calculated charge radii appear in good agreement with experimental ones, except 
in a few cases. One must stress in this respect that the charge radii obtained in the HFB 
approach do not include the influence of the long range correlations associated with collective 
oscillations of the mean-field, such as RPA ground state correlations or shape coexistence 
effects occurring in soft nuclei. As well known, these correlations beyond the mean field may 
lead to a significant increase of HFB r.m.s. radii. 

Fig. |l^ is the equivalent of Fig. ^ for neutron (upper ) and proton (lower part) radii. 
Experimental data on proton radii (crosses) have been taken from Ref. |ir3| , p!4| and the fits 
are given by formulae (|0|) and ( |6.2|) . 

Finally, as already pointed out in Section ||, the ratios Tp/ r„ are deformation- independent 
quantities that can be used to calculate the neutron radii of nuclei whose charged radii are 
known (see eq. ( p.l3| )). From formula ( p.5| ), one gets the following expression for neutron 
radii as functions of A, / = (A^ — Z)/A and (r^/j), the mean square charge radius: 



1.04(1 - 0.27/- 1.12/A) 



VII. CONCLUSIONS 

In this work, the results of self-consistent HFB calculations performed with the Gogny 
effective interaction for several isotopic and isotonic chains of spherical nuclei have been 
analysed and compared with those given by phenomenological expressions and by the rela- 
tivistic mean field approach. The more important conclusions that can be drawn from this 
analysis are the following: 

1. The average binding energy of nuclei calculated by applying the Strutinsky smoothing 
procedure to the single-particle level scheme obtained from the HFB self-consistent 
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mean field, eq. is in excellent agreement with the liquid-drop formula of My- 

ers and Swi§.tecki (M-S) [^]. The main difference occurs in the term correcting the 
charged liquid drop Coulomb energy for the diffuseness of the proton density. Also, 
the asymmetry coefficients in the volume and surface terms are found unequal, con- 
trary to the M-S formula, and the reduced radius parameter tq is slightly smaller (rg 
= 1.19 fm instead of tq = 1.205 fm). When all contributions are added, the differences 
between eq. ( ^.1|) and the M-S formula compensate, yielding nuclear binding energies 
in agreement within 3 MeV. 

2. The r.m.s. of the neutron and charge distributions agree in a quite satisfactory way 
with experimental data. From the set of nuclei studied in this work, systematics of 
the neutron, proton, charge and mass r.m.s. radii are obtained for spherical nuclei 
in the form of parametrizations depending on the mass number A and asymmetry 



/ = (A^ — Z)/A (see eqs. ( |6.1| )-( p7^ )). These parametrizations are of the same form as 
those first proposed in Ref. 0. The values of the parameters found in the present work 
are consistent with those derived from experimental data |^ and from the microscopic 
calculations performed in the framework of the relativistic mean field approach |Q. 

3. Systematics of the ratio rp/r„ of the proton r.m.s. to neutron r.m.s. radius are of 
special interest because they are, to a very good approximation, independent of the 
nuclear deformation and, in addition, they can be used to determine the value of the 
neutron radius of any nucleus, either spherical or deformed, when its charge radius 
is known (see eq. ( p.l3|) ). The parametrization of rp/r„ obtained from the set of 
nuclei considered in the present study (eq. (|6.5| )) appears in excellent agreement with 
previous ones based either on experimental data ^ or on the relativistic mean field 
approach 0. This parametrization can therefore be considered as robust enough to be 
used for predictions concerning nuclei for which experimental data not yet available, 
such as nuclei close to the neutron drip line or superheavy elements. 
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FIG. 1. {N,Z) plot of the spherical nuclei analysed in the present work. The circles indicate 
the /3-stable nuclei and the crosses a few representative isotope and isotone chains. 
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FIG. 2. The dependence on A of the parameter hcomin of the No=l4: HO bases employed in Ca, 
Sn and Pb isotopes. The error bars indicate ranges of fico for which the variation of the HF energy 
does not exceed 0.1 MeV. 
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FIG. 15. The fitted liquid drop energy E'ld given by eq. (solid lines) for constant Z (left), 

constant (center), and /3-stable (right) nuclei versus their mass number A. The differences 
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FIG. 16. Upper row: total HFB r.m.s. radii (solid lines), and their fits with formula ( |6.3f) 
(dashed lines) in constant Z (left), constant N (center) and /3-stable (right) nuclei. Lower row: 
HFB charge radii (dashed lines) compared with experimental data ||l^ (crosses) and with the fit 
given by formula ( |6.4D . 
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FIG. 17. Neutron (upper row, solid lines) and proton (lower row, dashed lines) HFB r.m.s. radii 
of constant Z (left), constant N (center), and /3-stable (left) nuclei compared with experimental 



data [|13|,U] (crosses) and with formulas (|6.1| )- (|6.2|) (dashed lines). 
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